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ABSTRACT 

It is now known that minor mergers are capable of creating structure in the phase- 
space distribution of their host galaxy's disc. In order to search for such imprints in the 
Milky Way, we analyse the SEGUE F/G-dwarf and the Schuster et al. (2006) stellar 
samples. We find similar features in these two completely independent stellar samples, 
consistent with the predictions of a Milky Way minor-merger event. We next apply 
the same analyses to high-resolution, idealised N-body simulations of the interaction 
between the Sagittarius dwarf galaxy and the Milky Way. The energy distributions of 
stellar particle samples in small spatial regions in the host disc reveal strong variations 
of structure with position. We find good matches to the observations for models with 
a mass of Sagittarius' dark matter halo progenitor S 10 11 M Q . Thus, we show that 
this kind of analysis could be used to provide unprecedentedly tight constraints on 
Sagittarius' orbital parameters, as well as place a lower limit on its mass. 

Key words: Galaxy: disc, structure - galaxies: formation - galaxies: kinematics and 
dynamics - methods: analytical - methods: iV-body simulations 



1 INTRODUCTION 

It is now widely accepted that the gravitational pull ex- 
erted by a merging satellite galaxy can generate struc- 
ture in the stellar disc of its more-massive host. A strik- 
ing example of the outcome of these kind of interactions 
is the system co mposed of M51 a nd its companion galaxy, 
NGC 5195 (e.g. lOh et all 120081 : iDobbs et all l20ld ). Pre- 
vious studies have shown that prominent spiral structure 
is more commonly observed in galaxies located within 
groups or with companions (Kormendv &: Normanl Il979l : 



lElmegreen fc Elmegreen 1 11982 . 1 19831 ). These observations 
highlighted the importance of tidal interactions on the struc- 



ture of galactic discs. Although grand design spirals, such 
as the ones observed in M51, are related to massive com- 
panions (in this case, Mmsi/Mngc 5195 > 0.3), several the- 
oretical studies have shown that much less massive satel- 



(Quinn, Hcrnquist & Fullagar 


19931: iTutukov & Fedorova 


2006: Kazantzidis et al.l 12008 


: Villalobos & Helmil 2008: 


Younger et al.l 20081: Quillen et al. 2009: Minchev et al. 


2009 


: Bird. Kazantzidis & Weinberg! l201ll; iGomez et al. 


2012 


). Morphological features, 


such as spiral arms, warps, 



flares, central bars, and low surface brightness ringlike struc- 
tures in the outskirts of a disc can be the result of either 
recent or ongoing accretion events. 
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Most of the previously-mentioned morphologi- 
cal features have been observed in the disc of our 
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own Galaxy ijLevine et all 2006; Mo manv et al.l 2006; 
Cabrera-Lavers et all 120071 ; iPohl et al.l 120081 ; iNewberg et al.l 



20021 : lYannv et al.l l2003l )T and various formation scenarios 
have been explored for most of them. An interesting 
example is the Monoceros ring, which is a low-latitude 
stellar structure spanning about 180° in Galactic longi- 
tude at nearly constant Galactocentric distance. It has 
been proposed that this structure could be the remnant 



of a tidally - disrup t ed satellite galaxy I Helmi et all 120031: 
Ibata et all 120031 ; IConn et al] 120051 ; IPenarrubia et al. I 



20051 ). More recently, several authors have argued that 



this dynamically cold ring could be the result of a tidal 
perturbation from a sa t ellite galaxy on the Galactic disc 
dKazantzidis et al.l [20081 : 1 Younger etai1l200Sl ; iQuillen et al] 
2009). In particular, the Sagittarius dwarf galaxy (Sgr) 
has been proposed as a possible perturber of the Milky 
Way, causing the emergence of the Mono ceros ring as 
the f ar extension of a Galactic spiral arm l|Purcell et al.l 
2011) and also potentially as a gravitational influence 
that circularised the orbit of a putative dwarf galaxy that 
could have been torn apa rt to become that same ring 
l|Michel-Dansac et al.l l201ll ). Motivate d by recent studies 
on co s mological abundance matching ([Conrov fc Wechsler 
120091; iBehroozi. Conrov fc Wechsler I 120100 Purcell et al.l 
l|201ll. hereafter Pll) considered models of Sagittarius 
with total masses as large as ~ 10% of the Milky Way's 
mass. These models can successfully reproduce not only 
dynamically cold structures such as the Monoceros ring, 
but also other global morphological features observed in 
the Milky Way. However, its impact on the phase-space 
structure of the Solar neighbourhood has not yet been 
explored. 

As first shown bv lMinchev et all ((2009), the energy kick 
imparted by the gravitational potential of a satellite as it 
crosses the plane of the disc can strongly perturb the ve- 
locity field of disc stars located in local volumes such as 
the Solar neighbourhood. In stellar samples close to the Sun 
(i.e. distances ^ 0.2 kpc), these perturbations can be ob- 
served in the u-v plane as arc-like features traveling in the 
direction of positive v, where u and v are the radial and 
tange ntial velocity compo nents, respectively. In a follow-up 
study. [Gomez et al.l (|2012l ) showed that the space defined by 
the energy and angular momentum of stars is a better choice 
than velocity space, as substructure remains visible even in 
much larger local volumes. Furthermore, they also showed 
that satellites with masses as small as 10% of the mass of 
the host can leave imprints in the phase-space distribution 
of Solar neighbourhood-like volumes that could be identified 
as late as w 5 Gyr after the first satellite's pericentre pas- 
sage. Substructure associated with this mechanism, known 
as "ringing," is expected to be better defined in the Galactic 
thick disc. This is because these stars spend relatively little 
time near the Galactic plane, where perturbations from the 
Galactic bar, spiral structure, and giant molecular clouds 
are more vigorous. Moreover, the thick disc is composed of 
a very old population of stars, with ages a round 10 to 12 Gyr 
l|Schuster et. al.ll2006l ; iBensbv et al.ll2007l ). Thus, most of its 
stars must have been in place at the time the hypothetical 
merger occurred. 

An obvious question arises from our previous discussion: 
Could a satellite galaxy like Sagittarius have left imprints of 
its tidal interaction with the Milky Way's thick disc in the 



local stellar phase-space distribution? Previous attempts to 
describe features in the velocity field of the Solar neighbour- 
hood have been primarily concerned with the dynamical ef- 
fects induced by non-axysimmetric disc components , such as 
a cen t ral bar or s e lf-gravitating spiral arms (see, e.g., Dehnen 
2000 1; iFuxl 120011; IMinchev et all 120071; IMinchev fc Quillen 
20081; lAntoia et al.l|2009l; iMinchev et al.ll2010l ; lAntoia et al. ' 



20 ill ; Quill en et al.ll201ir ). These studies have focused their 



attention on very small local volumes, dominated mostly by 
populations of thin-disc stars. However, little attention has 
been paid to the identification and characterisation of sub- 
structure in local volumes of the Milky Way's thick disc that 
may have originated as a respo nse of the tidal inter action 
with a satellite galaxy (see, e.g.. IMinchev et al]|2009h . 

In this work we attempt to fill this gap by analysing a 
full 6D phase-space catalog of F/G-type dwarf stars from 
the Sloan Extension for Galactic Understanding and Explo- 
ration (SEGUE; Yanny et al. 2009). The great advantage 
of this catalog, as presented bv lLee et al. I (|2011bl . hereafter 
Lll), is that it contains accurate estimates of [a/Fe] and 
[Fe/H] ratios for a large fraction of its stars. This enables 
a chemical separation of the disc system into likely thin- 
and thick-disc populations, avoiding unwanted biases asso- 
ciated with methods based on stellar kinematics or spatial 
distributions. To complement this a nalysis, we also explore 
the distribution of disc stars from the lSchuster et. all (|2006l ) 
catalog. 

The paper's outline is as follows. Section 2 presents a 
brief review of the main properties of ringing. In Section 3 
we describe the stellar samples analysed in this work. The 
distribution of stars in en ergy and angular mom entum space 
of the F /G-dwarf and the lSchuster et . all (|2006l) catalogs are 
analysed in Section 4 and 5, respectively. In Section 6 we ap- 
ply the same analyses to high-resolution N-body simulations 
of the interaction between Sagittarius and the Milky Way. 
A summary and conclusions are presented in Section 7. 



2 RINGING IN GALACTIC DISCS 



IMinchev et~ai] (|2009l ). followed bv lGomez et all l|2012l . here- 
after G12), showed that relatively massive minor-merger 
events can generate substructure in the velocity field of disc 
stars located in the Solar neighbourhood. Merging satel- 
lite galaxies with total masses as small as 10% of the mass 
of the Milky Way may excite density waves whose signa- 
tures can be identified even 5 Gyr after the first satellite's 
pericentre passage. The space of E-L z is particularly well- 
suited to identify density waves, since the waves are sur- 
faces of constant energ y (see also e.g. iHelmi fc de Zeeuwl 
120001 ; iHelmi et al.lll999r ). Over time, these constant-energy 
features associated with the density waves become more 
closely spaced, and their number increases as a consequence 
of phase wrapping (see Section 3 of G12). In the top panel of 
Figure[T]we show, in E-L z space, the distribution of disc par- 
ticles located within a Solar neighbourhood-like sphere of 2.5 
kpc radius, obtained from one of the simulations analysed in 
G12. The sphere is centred at 8 kpc from the galactic centre 
and the distribution was obtained after 3.9 Gyr of evolu- 
tion, which corresponds to a time since total disruption of 
the satellite galaxy of 1.9 Gyr. In this simulation, the satel- 
lite galaxy has a total initial mass equal to 20% of the mass 
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Figure 1. Top panel: E — L z distribution of disc particles in- 
side a sphere of 2.5 kpc radius located at 8 kpc from the galac- 
tic centre, obtained from a simulation with a mass ratio of 
^hostV*^sat al = °- 2 - The total number of particles enclosed 
within this volume is N™° v ^ al = 976. The snapshot corresponds 
to a time of 3.9 Gyr after the satellite was launched, or to 1.9 Gyr 
after the satellite fully merged with the host galaxy. The different 
colours (contours) indicate different number of particles. Note the 
large number of lumps with nearly constant energies associated 
with density waves crossing the sphere. Bottom panel: Kernel his- 
togram (KH) of energies, E, obtained from the same distribution 
of particles. Each peak in the kernel histogram corresponds to a 
density wave observed in the top panel. 



of the host. Note that only N™?* 1 = 976 are located within 
this 2.5 kpc sphere. It is clear that disc particles are dis- 
tributed in different lumps of nearly constant energy, each of 
them associated with a different density-wave crossing. The 
bottom panel of Figure [1] shows a Kernel Histogram (KH) of 
energies, E. To compute the kernel histogram, a Gaussian 
kernel with bandwidth a = 0.011 AE was assigned to each 
particle. Here, AE — 4 x 10 4 km 2 s -2 is approximately the 
total extent of the satellite in energy space within this vol- 
ume. The kernel histogram is the sum of all of the Gaussian 
kernels. Since density waves are surfaces of nearly constant 
energy, they are observed in the KH as well-defined peaks. 
Note that larger values of a tend to oversmooth the KH, 
erasing some of the peaks associated with substructure ob- 
served in E — L z space. Conversely, smaller values produce 



noisier KHs. In what follows we keep the values of a and 
AE fixed at these values. 



3 THE STELLAR SAMPLES 

3.1 The SEGUE F/G-Dwarf sample 

The SEGUE F/G-dwarf sample analysed in this work was 
culled from ~ 70, 000 F and G type stars with avail- 
able low-resolut ion spectroscopy (R ~ 2000), as provided 
in SDSS DR8 jAihara et all l201lh . Among them, about 
63, 000 stars were targeted as G-dwarfs; hence, our sample 
is dominated by the G-dwarf candidates used by Lll. The 
G-dwarf candidates were obtained by selecting stars with 
colours and magnitudes in the range 0.48 < (g — r)o < 0.55 
and ro < 20.2, respectively, while the rest of the sample 
covers 0.2 < (g - r) < 0.48, and g < 20.2. Thanks 
to this simple selection function, the stellar sample is ex- 
pected to be completely unbiased with respect to kinemat- 
ics, and only slightly b iased towards metal poor stars (see 
ISchlesinger et alj|2009l 1. 



3.1.1 Local sample 

In order to obtain a local sample of stars with accurate mea- 
surements of their 6D phase-space coordinates, as well as 
metallicities and [a/Fe] ratios, a series of cuts were applied 
to the complete catalog. Here we provide a brief summary 
of the criteria used for the sample selection, but we refer the 
interested reader to Lll for a more detailed description of 
this procedure as we follow their prescription to select the 
final sample. 

Stellar atmospheric parameters, such as effective tem- 
perature, T a a, surface gravity, log g, and metallicity, 
[Fe/H], were deter mined using the S E GUE Stellar Parame- 
ter Pipeline, SSPP (TLee et al.ll2008al lbl: lAllende Prieto et all 
120081 : ISmolinski et al. 2011). In order to obtain high-q uality 
estimates of both [Fe/H] and [a/Fe] (|Lee et al.ll2011aT ) . only 
stars with spectra of signal-to-noise ratios (S/N) greater 
than 30 A were considered. Importantly, this also ensures 
that errors in the estimated radial velocities are smaller 
than 5 km s -1 . P roper motions were obtained following 
iMunn et al.l ll2004h. after correctin g for the systematic er- 
ror described in Munn et alj (|200&t ). Distances to individual 
stars were estimat ed using calibrated set of stellar i sochrones 
An et ail l2009rj ). following the prescription of lAn et all 



2009aj). To minimise possible distance bias from stellar 



age effects near the main sequence turnoff, only stars with 
log g ^ 4.2 were considered. In addition, to minimise the 
errors on the estimates of phase-space coordinates, we only 
consider stars with distances from the Sun 0.4 < d ^ 2 kpc. 
Note that the inner limit is imposed by the data itself, due 
to saturation of the SDSS photometric scans for g < 14.5. 

To compute Galactocentric positions and velocities of 
our local stellar sample, we assume the Sun to be located 
at T?© = 8 kpc, and that it has a peculiar velocity with re- 
spect to the Local Stand ard of Rest (LSR) (U , V, W)q = 
(11.1, 12.2, 7.3) km s _1 (|Schonrich et al.ll2010D . We further 
assume a velocity of the LSR with respect to the Galactic 
centre of Vlsr = 220 km s . Finally, we discard stars with 
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[Fe/H] ^ — 1.2 and V$ < km s _1 in order to avoid con- 
tamination (as much as possible) from the stellar halo, as 
well as t he proposed metal-w eak thick-disc component of the 
Galaxy (|Carollo et al.ll201Gft . 



3.1.2 The thin- and thick-disc populations 

As mentioned in the Introduction, signatures of ringing are 
expected to be stronger in the thick disc than the thin 
disc because its stars are less affected by perturbations 
from the Galactic bar, self-gravitating spiral arms, and gi- 
ant molecular clouds. More importantly, the thick disc is 
mainly composed of a very old stellar population, with age s 
of - 10-12 Gyr (jSchuster et. al.ll200a Bensbv et al.ll2007h . 
Thus, most of its stars must have been in place as recently as 
5 Gyr ago - a required condition if, for example, we seek to 
identify signatures of a merger event that may have started 
approximately at that time. 

Following Lll, we split our local sample into likely thin- 
and thick-disc populations on the basis of their stellar chem- 
ical abundances, i.e., [a/Fe] and [Fe/H]. We assign stars to 
the thin-disc (low- [a/Fe]) and thick-disc (high- [a/Fe]) com- 
ponents according to the following scheme: 



For stars with [Fe/H] ^ -0.: 

• thin disc, if [a/Fe] < - 

• thick disc, if [a/Fe] > - 



0.08 • 
-0.08 



[Fe/H] +0.15 
■ [Fe/H] +0.25 



For stars with [Fe/H] < -0.8 

• thin disc, if [a/Fe] < +0.214 

• thick disc, if [a/Fe] > +0.314 

While any old stellar population with high velocity dis- 
persion would be suitable for our investigation, by assigning 
stars to the different components of the disc according to 
[a/Fe] we avoid introducing unwanted biases associated with 
methods based on stellar kinematics or spatial distributions. 
Note the gap of 0.1 dex left in the [a/Fe] cuts between the 
thin disc and thick disc, which is chosen to avoid misclassifi- 
cation of stars. For a detailed description of the motivation 
behind this scheme, we refer the reader to section 3.1 of Lll. 



3.2 The Schuster et al. (2006) sample 

The ISchuster et. al.1 (|2006l . hereafter SCH06) catalog com- 
prises a total of 1533 high-velocity and metal-poor stars. 
Positions and proper motions for these stars were pr i- 
marily d erived from Hi pparcos (|Perrvman et al.l Il997h . 
Tvcho-2 toig et al.ll2000l ). and the revised NLTT catalog 
l|Salim fc Gouldl 120031 ). whereas radial velocities were ob- 
tained from a number of different sources available in the 
literature (see SCH06, and references therein). As in Sec- 
tion 13.11 wc compute Galactocentric velocities by assuming 
R e = 8 kpc, (U, V, W) e = (11.1, 12.2, 7.3) km s' 1 and 
VisH = 220 km s" 1 . 

Metallicities estimated via a photometric calibration 
based on uvby-/3 photometry are also provided for the entire 
sample. As shown by SCH06, it is possible to discriminate 
stars from different Galactic components by means of the X 
parameter, where X is a simple linear combination of the 
rotational velocity and metallicity of the stars. Following 
SCH06, we discard from our sample all stars with X ^ 0, 



Table 1. Parameters of the present day Milky Way-like potential used 
in this work. Masses are in Mq and distances are in kpc. 



Disc 


Bulge 


Halo 


M d = 4 x 10 10 


M b = 8 x 10 9 


M vir = 1 x 10 12 


r a = 6.5 


r c = 0.7 


r s = 21.5 


r b = 0.26 




c = 12 



in order to avoid potential contamination from the Galactic 
halcfl 

Note that estimates of [a/Fe] are not provided in the 
SCH06 catalog. Thus, a subdivision between likely thin- 
and thick-disc stars cannot be performed following the 
sc heme described in Section 13.1.21 Nevertheless, as shown 
bv lNissen fc Schusterl (|l99ll ) and SCH06, the thin-disc pop- 
ulation in this sample belongs to a very old disc component, 
with the bulk of the population older than 10 Gyr, and the 
remaining stars with ages between 10 and 4 Gyr. We thus 
ensure that most of these stars were already part of the disc 
~ 5 Gyr ago. It is important to note that all disc stars in this 
sample are located at a distance from the Sun d < 0.2 kpc. 
Thus, this sample probes the innermost regions of the So- 
lar neighbourhood, complementing the more distant SDSS 
F/G-dwarf sample. 



4 THE SEGUE F/G-DWARF SAMPLE IN 
E - L z SPACE 

We now analyse the E — L z distribution of the stars present 
in our local SEGUE F/G-dwarf sample. In order to obtain 
an estimate of the orbital energy of these stars we first need 
to assume a shape for the Galactic potential. We adopt a 
model consisting of three dif ferent components. Thes e are: 
(1) a Miyamoto-Nagai disc (|Mivamoto fc Nagaill"l975h 



<&disc = 



GMd 



(2) a Hernquist bulge (|Hernqu ist 1990), 



$bul. 



GMb 
r + r c 



(1) 



(2) 



and (3) a NFW dark matter halo ()Navarro. Frenk fc White! 
1996) 



GM vir 



r [log(l + c) - c/(l + c)] 



log 



(3) 



Table [T] summarises the numerical values of the pa- 
rameters used for our Galactic potential model. These 
parameters are similar to those used by other au- 
thors (e.g.. iBullock fc Johnston! 120051 : iGomez et al] l20ld : 
iPefiarrubia et al. 1 120101) a nd fit the Milky Way rotation 
curve ilKlypin et all 12002ft . Note that, as described by 



IGomez fc Helmil (|2QloT and references therein), the results 
are not strongly dependent on the particular choice of the 
potential, since substructure in integrals of motion space are 

1 Note that only a fraction smaller than 3% of our final F/G- 
dwarf thick-disc subsample has X ^ 0. 
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Figure 2. Top panel: E — L z distribution of likely thick-disc stars 
located within 2 kpc of the Sun, obtained from our F/G-dwarf 
sample. The different colours (contours) indicate different num- 
ber of stars. Only stars with total errors in velocity V eTT < 20 
km s" 1 are considered. The subsample contains iV'£j£~ = 3141 
high- [cr/Fe] stars. Note the presence of a series of high-density 
features, especially at large values of energies. Bottom panel: ker- 
nel histogram of E, obtained from both the thick-disc (red solid 
line) and the thin-disc (blue dashed line) stellar subsamples. The 
thin-disc subsample contains JV^JjJ^ = 4903 low-[a/Fe] stars. The 
shaded areas indicate errors associated with poor sampling of 
the underlying distributions (see text). The high-density features 
observed in the top panel correspond to peaks in the kernel his- 
togram. Note that the thin-disc subsample exhibits a smoother 
distribution relative to the thick disc. 



robust to small differences in the mass distribution because 
we always focus on small volumes in space. Hence, small 
changes in the potential essentially act as a zero-point off- 
set, affecting all the stars present in this volume in the same 
way. 



figure we have only considered stars with d ^ 2 kpc from 
the Sun and a total velocity error V err < 20 km s . As a 
result, and because of the different cuts applied to the com- 
plete catalog (see Section f3.1|l . we are only left with a total 
of AS? = 3141 likely thick-disc stars. A series of high- 
density features can be observed in this panel, especially at 
high values of E and L z . These features become more evi- 
dent in the kernel histogram of energies, shown as a red solid 
line in the bottom panel of the same figure. Comparison with 
Figure [1] clearly indicates that, if these features are to be as- 
sociated with ringing, they have a much smaller amplitude 
than might be expected. However, due to the relatively large 
threshold in total velocity errors considered (i.e., V elI < 20 
km s -1 ), these peaks may have been smoothed out; we ex- 
plore this issue further in Section 4.2. On the other hand, 
due to the relatively low number of stars, some of the peaks 
may have arisen due to noise associated with poor sampling 
of the underlying distribution. To investigate this, we boot- 
strapped our thick-disc stellar subsample 2000 times and 
computed a kernel histogram for each one of the realisa- 
tions. At each value of E, we sort the results of the KHs in 
ascending order. The shaded areas in the bottom panel of 
Figure [2] show the 25-75 percentiles (light shadow) and the 
5-95 percentiles (dark shadow) of the distribution of values 
obtained from the bootstrapped KHs. The results of this 
procedure indicate that only a few of these peaks have am- 
plitudes with values above the noise level. However, in the 
following sections we show that, after accounting for the ef- 
fects of measurement errors, and through exploration of the 
SCH06 catalog, some of the remaining peaks may indeed be 
real. 

The bottom panel of Figure also shows (blue dashed 
line) the kernel histogram on E of the stars associated with 
the low-a, likely thin-disc stars. In this case, the total num- 
ber of stars in the subsample is N^rs — 4903. It is very 
interesting to observe that the distributions associated with 
the thin disc and thick disc are shifted with respect to one 
another. The thin-disc subsample comprises a population of 
stars that, at least locally, is less tightly bound than that of 
the thick disc. Note that, within volumes as small as the ones 
considered here, the Galactic potential can be regarded as 
constant. Thus, for these stars, the binding energy is mainly 
a function of their velocities. Since the thick disc is com- 
posed of a population of stars with a lag in V$ of ~ 30 km 
s _1 relative to the thin disc (as shown in, e.g., LIT), it is 
not surprising to find that these stars have a larger binding 
energjQ. It is also noticeable that the thick-disc stars exhibit 
a much broader distribution in energy than their thin-disc 
counterparts. This is a manifestation of the larger velocity 
dispersion (in all directions) associated with the thick-disc 
component. Finally, note that the thin disc exhibits a much 
smoother distribution when compared to that of the thick 
disc. One has to bear in mind that the thin-disc subsample is 
also ~ 35% larger than the thick-disc subsample. Thus, noise 
associated with poor sampling is reduced. This is indicated 
by the tighter blue and light-blue shaded regions, computed 
as discussed above. Nevertheless, in Section 14.21 we show 



4.1 Thin vs. Thick disc 

The top panel of Figure [2] shows the distribution of a- 
enhanced, likely thick-disc stars in E ~ L z space. For this 



2 Note that the mean radial and vertical velocity components for 
both populations are ~ km s . 
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Figure 4. Kernel histograms of normalised energies, obtained from the distribution of particles shown in Figure [T] after convolving 
with measurement errors. The red solid line in each panel shows the averaged kernel histogram obtained after 400 realisations of the 
error convolution. The maximum error in total velocity considered is indicated in the top right corner of each panel. The shaded regions 
indicate la deviations from the mean. For comparison we show, with a dashed black line, the kernel histogram obtained from the same 
distribution of particles before convolution. When errors in total velocity as large as 20 km/s are considered, most of the peaks in the 
histogram are erased (left panel). These peaks are recovered when we reduce the magnitude of the maximum error to 11 km/s (right 
panel). 
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Figure 3. Distribution of total errors in positions (top panel) 
and velocities (bottom panel), extracted from our likely thick- 
disc SDSS F/G-dwarf subsample. 



that this result holds even when considering smaller stellar 
subsamples with more accurately determined velocities. 



4.2 Characterising the effects of errors on 
phase-space coordinates 

As noted above, the measurement errors in the phase-space 
coordinates of our local F/G-dwarf stellar sample may be 
large enough to either smooth out or erase some of the sig- 
natures of ringing in E — L z space. To characterise their 
effect on the observed features in the kernel histograms, we 
have extracted the distributions of errors from our thick- 
disc subsample, and convolved them with the position and 
velocities of the particles shown in Figure [T] Estimates 
of the errors were obtained from the SSPP (see Lll and 
references therein); their distributions are shown in Fig- 
ure [3] The left panel of Figure [4] shows the results of this 
convolution when errors in velocities as large as 20 km -1 
are considered. The red solid line shows the averaged ker- 
nel histogram of normalised energies, -Enorm, obtained after 
400 different realisations of the error convolution, whereas 
the shaded region indicates its standard deviation. Here, 
£ Eorm = (Ei — mm(E))/AE, with AE defined as in Sec- 
tion [2] It is important to remark that, to compute the en- 
ergies of the iV-body particles after convolution with errors, 
we have used the same analytic potential that was used for 
the F/G-dwarf sample. Note that this simulation was de- 
signed to emulate the result of a minor merger experienced 
by a Milky Way-like galaxy at z — 1. Thus, the properties 
of this potential were scaled to those expected at z = 1, as 
described in Section 2 of G12. Interestingly, due to the large 
magnitude of the errors considered, most of the structure 
associated with ringing has been erased. Furthermore, some 
spurious peaks can be observed. For comparison we show 
with a black dashed line the kernel histogram obtained from 
this set of particles before error convolution. The middle 
panel of this figure shows the results of the same procedure, 
now after convolving with errors in velocity Verr < 16 km 
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s _1 . Note that these errors are still large enough to erase 
most of the signatures of ringing. Only when considering er- 
rors in velocity Vi rr < 11 km s _1 (right panel) do we start 
recovering the previously-observed peaks. 

This exercise clearly indicates that, in order to detect 
signatures of ringing in the Milky Way disc, very accurate 
6D phase-space catalogs must be explored. A highly accurate 
subset of stars can be obtained from our local F/G-dwarf 
sample, at the expense of significantly reducing the total 
number of stars analysed. The top panel of Figure [S] shows 
the kernel histograms of the thick disc (red solid line) and 
the thin disc (blue dashed line) obtained when only stars 
with total velocity errors V crr < 11 km s _1 are considered. 
As a result of this more restrictive cut, we are left with 
iVS = 968 thick-disc and N^ s = 1132 thin-disc stars, 
respectively. Interestingly, not only do we recover most of 
the features already observed in Figure [2] but also structure 
here is more sharply defined (as expected from our previous 
discussion). Note, however, that errors associated with poor 
sampling of the underlying distribution are larger due to the 
smaller number of stars. This is evident by the wider areas 
covered by the shaded regions, obtained from the bootstrap 
analysis of our subsample (see Section 14. ip . Nevertheless, 
it is important to note that the number of thick-disc stars 
in this more accurate subsample, iVg t ^f, is approximately 
equal to the number of particles analysed in the iV-body 
model, N££t el , shown in Figures [T] and 0] 

Alternately, it is possible to obtain a more accurate sub- 
sample simply by analysing a much more local subset of 
stars. In the bottom panel of Figure [5] we show the kernel 
histogram of E obtained from a sample of stars with he- 
liocentric distances projected onto the plane R ^ 0.5 kpc. 
Note that by reducing the radius of our Solar neighbourhood 
sphere we are not only constraining errors in the velocities of 
the stars, but also in their positions. After the cut we obtain 
a thick-disc subsample of N^^s = 2441 stars, whereas the 
thin-disc subsample contains A^ t h a '" s = 3021 stars. Compari- 
son with the top panel of Figure [5] and Figure [2] shows again 
that many previously-observed features are preserved. Fur- 
thermore, some of the features have become even more sig- 
nificant, in particular, the peak located at E » —0.97 x 10 
km 2 s -2 . It is also interesting to note that, as previously 
discussed in Section T4. II the thin-disc population exhibits a 
smoother energy distribution relative to that of the thick- 
disc population in both of these panels. Notice that now 
both subsamples have comparable numbers of stars. 



5 THE SCHUSTER ET AL. (2006) SAMPLE IN 
E-L z SPACE 

The previously-observed peaks in the kernel histograms of 
our SEGUE F/G-dwarf sample might be an indication of 
ringing in our own Galactic thick disc. However, due to the 
relatively small number of stars with the necessary small 
measurement errors, uncertainties associated with poor sam- 
pling of the underlying distribution are (for some features) 
still large. It is therefore important to look for additional 
evidences of ringing in other stellar catalogs already avail- 
able in the literature. The SCH06 is particularly well-suited 
for this analysis, as it comprises a set of old and metal-poor 
disc stars (see Section [3.2(1 . Furthermore, the disc stars in 
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Figure 5. Top panel: Kernel histograms of energies obtained from 
the SEGUE F/G-dwarf thick disc (red solid line) and the thin 
disc (blue dashed line), when errors in total velocity smaller than 
11 km/s are considered. The shaded regions indicate errors as- 
sociated with poor sampling of the underlying distribution (see 
Section 14.11 1 . Note that most of the previously-observed features 
in Figure f2] are more sharply defined, as expected from this more 
accurate subsample. Note as well that the thin disc presents a 
smoother distribution in -Enorm with respect to the thick disc. 
Bottom panel: As above, for stars located within R ^ 0.5 kpc 
and total velocity errors smaller than 20 km/s. 

this sample are all located within d < 0.2 kpc, thus they 
probe a region of the Solar neighbourhood that is different 
than that explored by the more distant SEGUE F/G-dwarf 
sample. 

The black dashed line on the top panel of Figure[6]shows 
the kernel histogram obtained from a subsample of likely 
disc stars from the SCH06 catalog. Following our previous 
discussion, we first consider stars with total errors in veloc- 
ities Verr ^11 km s . As a consequence, we are left with 
a total of 603 old disc stars. Interestingly, a large number 
of peaks can be observed, most of them located at values 



© RAS, MNRAS 000,[TJfl4] 



0.75 



o 

X a 0.5 

2 



0.25 



A. Gomez et al. 




V < 1 1 km/s ' 

err 


M 


1 ^ 
I 'l 


1 i ii 
1 i 


1 1 I 
1 ' I 
l> I 

V 1 

V i 

\ i 


TV /v' 

A / / 
/ \ / ' 
1 \l ' 
1 ' 
/ 

y " 


V\ V ' 

\ l 

1 • /i 

\ >' 1 


-1.2 -1.1 


-1 -0.9 




£ [10 5 km 2 ^ 



Figure 6. Top panel: Kernel histogram of E obtained from the 
SCH06 sample of disc stars with total velocity error V eTT < 11 
km/s, shown with a black dashed line. The sample contains a 
total of 603 stars. Note the large number of peaks located at 
very similar values of E as in the kernel histogram of the F/G- 
dwarf subsample, indicated with a red solid line. The slight shift 
between these two distributions can be accounted for by the dif- 
ferent volumes probed by the samples (see text), as well as poor 
sampling of the underlying energy distribution. Bottom panel: 
As in the top panel, for a sample of disc stars with total veloc- 
ity error V er r < 20 km/s. This larger sample contains 813 stars. 
The coloured areas indicate the approximate energy range of the 
Arcturus stream (yellow), the Hercules stream (green), and the 
most significant peak identified in the SEGUE F/G-dwarf sample 
(light blue). 



of E consistent with those observed in the kernel histogram 
of F/G-dwarf subsample, shown with a red solid line. As 
we explain in Section T6.2I due to the strong radial and az- 
imuthal dependance on the location and amplitude of the 
peaks, together with the poor sampling of the underlying 
energy distribution, it is not surprising to find features in 



these kernel histograms slightly shifted with respect to one 
another. 

The bottom panel of Figure [6] shows the kernel his- 
togram from the SCH06 catalog, obtained after considering 
stars with total velocity errors of Krr < 20 km s" 1 . This 
larger sample contains 813 stars. Again, note that many of 
the features observed in the KHs of both the SCH06 and the 
F/G-dwarf samples are located at very similar values of E. 
The colour-coded areas in this panel indicate the approxi- 
mate energy range in which the Arcturus and the Hercules 
moving groups are expected to be found. To compute these 
energy ranges, for simplicity, we have assumed for each mov- 
ing group: 

• a charact er istic V vel o city c o mponent, Vls r (see , 
e.g.. iDehnenl 1 19981: iFuxl |200~D : lAntoia et al.l 120081: 
[ Navarro. Helmi fe Freeman! I2004E iKlement. Fuchs fc Rbd 
|2008| ). Since this characteristic V velocity is subject to 
relatively large uncertainties, to compute each energy range 
we assumed V = Vlsr ± 10 km s _1 

• a characteristic U velocity component equal to zero 

• a Galactocentric distance R = 8 kpc 

The bottom panel shows that both moving groups could 
be associated with peaks in our KH. Note that, as previ- 
ously explained in this Section, it is not surprising to find 
slight shifts in the location of the peaks between the two 
analysed samples due to the different volumes probed. M09 
considered the possibility that the Arcturus moving group 
is the result of a perturbation in the Galactic disc induced 
by a minor merger event. They concluded that a minor 
merger that took place ~ 1.9 Gyr ago could explain not 
only the presence of the Arcturus moving group, but also 
other pr eviously observed features in the local velocity field 
(see e.g. lArifvanto fc Fuchsll2006h . It is important to note 
that streams with velocities larger than Vlsr ~ ±50 km s _1 
are less likely to be related to the Galactic bar. Interestingly, 
the short timescale associated with this event is consistent 
with the estimated age o f the Galactic bar, as measured by 
ICole fc Weinberg I (|2002r i (< 3 Gyr ago). This suggests that 
the same event that could have caused the formation of the 
Galactic bar could have also left the stellar disc unrelaxed, 
thus giving rise to the observed high- velocity stream. On 
the other hand, many studies have successfully described 
the presence of the Hercules moving group as a result of 
a resonant interactio n with the bar (see e.g. lDehnenll2000l; 
iMinchev et al.l I2007T) . However, as shown by Antoia et al.l 
(2009. 1201 lh . the effect that bars and spiral density waves 
have on the velocity field of galactic discs is highly degener- 
Therefore, the scenario in which the Hercules moving 
group was formed as the response of the Galactic disc to a 
minor merger event cannot be ruled out. 



6 THE SAGITTARIUS DWARF GALAXY AS A 
POSSIBLE PERTURBER 

The impact that the Sagittarius dwarf galaxy may have 
had on the Galactic disc has been previously considered 
by several authors. Recently, using high-resolution Af-body 



3 Note that both mechanisms could be acting simultaneously. 
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Table 2. Initial properties of the two iV-body simulations analysed 



in Section 6. 


Host 


DM halo 






-Wpart 


= 2.65 x 10 7 


Virial mass 
Scale radius 


1 X 10 12 
14.4 




[Mq] 
[kpc] 




Stellar disc 






-Wpart 


= 3 x 10 6 


Mass 

Scale length 
Scale height 


3.59 x 10 10 

2.84 

0.43 




[M e ] 
[kpc] 
[kpc] 




Stellar bulge 






-Wpart 


= 5 x 10 5 


Mass 

Effective radius 
Sersic index 


9.52 x 10 9 

0.56 

1.28 




[M ] 
[kpc] 




Satellites 


DM halo 


Light 


Heavy 


-Wpart 


= 1.8 x 10 6 


Virial mass 
Scale radius 


0.32 X 10 11 
4.9 


1 X 10 11 
6.5 


[M ] 
[kpc] 




Stellar spheroid 






-Wpart 


= 5 x 10 4 


Core radius 
Tidal radius 
Central vel. disp. 


1.5 

4 
23 


1.5 

4 

30 


[kpc] 
[kpc] 
[km s 


- 1 ] 



simulations, Pll showed that its gravitational interaction 
may play a significant role on shaping the morphology 
of the Galactic disc, by inducing the formation of rings, 
influencing the Galactic bar, and flaring the outer disc. 
Its interaction c ould also explain the formation of the 
Monoceros rings llNewberg et al.ll2002l: iQuillen et al.l [20091 : 



iMichel-Dansac et all 120111 : iPurcell et all 120111 ). However, 
until now, its influence on the phase-space distribution of 
disc stars located within the neighbourhood of the Sun has 
not been explored. In this section we analyse the simulations 
presented by Pll to characterise the effect of Sagittarius 
within this volume. 



6.1 The simulations 

Two simulations with different models for the Sagittarius 
dwarf galaxy progenitor were performed. In both cases, the 
primary galaxy includes a NFW dark matter (DM) halo, 
an exponential stellar disc and a central bulge following 
a Sersic profile. The Light Sgr and Heavy Sgr progenitors 
were self-consistently initialised with a NFW DM hal o and 
a sep arate stellar component following a King profile (|Kind 
1966). Tabled summarises the numerical values of the pa- 
rameters used for t hese models. Following previous wo rk on 
the Sgr interaction (|Keselman. Nusser fc Pee bles 2009]), the 
satellites were launched at 80 kpc from the galactic centre 
in the plane of the Milky Way, traveling vertically at 80 km 
s _1 toward the north galactic pole. The simulations reach 
a present day configuration after approximately 2.7 and 2.1 
Gyr of evolution, respectively. Mass loss that would have oc- 
curred between virial radius infall and this "initial" location 
is accounted for by truncating the progenitor DM halo mass 



profile at the instantaneous Jacobi tidal radius, r t = 23.2 
and 30.6 kpc, respectively. This leaves a total bound mass 
that is a factor of ~ 3 smaller than their effective virial mass 
originally assigned. Each of the model Sgr progenitors expe- 
riences two disc crossings, approaching a third at the present 
day. The satellite first crosses the disc at a galactocentric dis- 
tance of ~ 20 kpc approximately ~ 1.75 Gyr ago, produc- 
ing the most significant perturbation. The progenitor loses 
roughly 75% of its dark matter mass (but little stellar ma- 
terial) during this time. In the bottom panel of Figure [7] we 
show the time evolution of the satellites' galactocentric dis- 
tance. The initial disc in both Sgr-infall models is completely 
smooth at t = Gyr, and only develops a mild bar when it 
is evolved in isolation for a few Gyr. When compared to the 
isolated run, the Light Sgr model results in a stronger bar 
(see Pll). Conversely, as a result of enhanced central disc 
heating, the more massive satellite suppresses bar formation. 
The endstate bar orientation, <j)^™ s = 15° - 20°, is in both 
cases consistent with estimates of the long bar's orientation 
observed in the center of t he Milky Way, = 15° - 30° 

(Biss antz fc Gerhardll2002T ). The simulations used the paral- 
lel A^-body tree code ChaNGa (|Gioachin et aUl2007h . with 
a gravitational softening length of one parsecQ, and followed 
the evolution of 30 million particles with masses in the range 
1.1 - 1.9 x 1O 4 M . 



6.2 The phase-space distribution of Solar 
neighbourhood-like volumes 

We now characterise the phase-space distribution of parti- 
cles within Solar neighbourhood-like volumes of 1 kpc radii. 
In these simulations, according to the location of the Sagit- 
tarius remnant and its tidal debris, the Sun should be lo- 
cated at approximately (0, —8, 0) kpc from the galactic cen- 
tre. In order to examine the dependencies of the particle's 
phase-space distribution with azimuthal angle, and uncer- 
tainties on the location of the "Sun" associated with the 
model, we have placed 5 spheres at different azimuthal an- 
gles, covering a total of fa 58° around the "Sun." The top 
panel of Figure [7| shows contour plots of the final distribu- 
tion of particles in the X - Y plane obtained from the Heavy 
Sgr simulation. The white dots trace the regions where Sgr 
crosses the plane of the disc and its current projection on 
the X- Y plane. The locations of our Solar neighbourhood 
spheres are indicated with solar symbols. As shown by Pll, 
the gravitational interaction between the Sagittarius-like 
satellite and the galactic disc is sufficiently strong to in- 
duce the formation of transient spiral arms or rings. These 
arms can be clearly observed outside a radius of 15 kpc, 
but not so clearly at the Solar radius. In the top panels of 
Figure [S] we show the kernel histograms of E obtained from 
the distribution of particles located within these volumes. 
Let us recall that throughout this work we have kept fixed 
both the value of AE and the value of the bandwidth of 
the Gaussian kernel, a. The black solid lines show the re- 
sults obtained when the particle energies are computed using 
the analytic potential introduced in Section [J] whereas for 
the red dashed lines energies were computed directly from 



The force softening length is set to roughly 0.1 of the mean 
disc interparticle spacing. 
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Figure 8. Top panel: Kernel histograms of energies obtained from the Heavy Sgr simulation. Each panel shows the kernel histogram 
(KH) obtained from the distribution of particles within a Solar neighbourhood-like sphere located at a different galactocentric azimuthal 
angle, 8. The black solid lines show KHs obtained after computing particle's energy using the analytic potential described in Section [4] 
The red dashed lines show the KHs obtained after computing particle's energies from the JV-body potential. The shaded regions indicate 
errors associated with poor sampling of the underlying distribution (see Section |4Tj . The angle 8 is measured clockwise, starting from 
the negative Y axis in Fig. [7] The azimuthal location of the spheres are indicated in the top left corner of each panel. Note that all the 
KHs exhibit multiple peaks associated with density waves crossing the Solar neighbourhood-like volumes. Furthermore, the locations of 
these peaks strongly depend on the azimuthal location of the sphere. Note as well the good agreement between the second panel (8 Ri 18° 
) and the KHs obtained for our F/G-dwarf thick-disc subsample shown in Figures [2] and [5] Middle panels: As above for the Light Sgr 
simulation. Note that, as expected for a less-massive satellite, the peaks previously observed on the KHs of the Heavy Sgr model are 
either very mild or not present at all. Bottom panels: As above, now for spheres located at a different galactocentric radius, R. The 
spheres are located at an azimuthal angles = 0°. The galactocentric distance of the spheres are indicated in the top left corner of each 
panel. The black solid and blue dashed lines show the KHs obtained from the from the Heavy Sgr and light Sgr simulations, respectively. 
For the Heavy Sgr simulation, the locations and amplitude of the peaks in the KHs strongly depend on the radial position of the spheres. 
Note as well the similarities between the second panel (R = 7.5 kpc) and the KHs obtained for our F/G-dwarf thick-disc subsample. For 
the Light Sgr simulation, mild signatures of density waves can be observed only in the outer spheres, where the perturbation from the 
satellite galaxy is stronger. 



the JV-body potential. The very good agreement between 
the two KHs reflects the weak dependence of the amount 
and distribution of features with the particular choice of the 
galactic potential. Note that, for comparison, the red dashed 
lines have been slightly shifted in energy, so that both KHs 
overlap in the same energy range. On average, the number 
of particles contained in each sphere is ~ 10, 000. It is in- 
teresting to observe that all the analysed volumes exhibit 
peaks in their KHs; these are features associated with den- 
sity waves. We see strong dependence of the distribution and 
amount of peaks with azimuthal angle, 9: as 9 is decreased, 
power from the highest amplitude peak (E w —1.01 x 10 5 
km 2 s -2 at 8 rss 29°), or density wave, is transferred towards 
a secondary peak at higher energy (E » —0.96 x 10 s km 2 



s -2 at 9 ~ —29°). In addition, the location of each peak is 
shifted towards lower values of E indicating that, on aver- 
age, these peaks are populated by more inner particles. This 
is the expected behaviour from a spiral pattern traveling in 
this angular direction (see also Section 4.2 of G12). Note 
the similarities between the KH in the second top panel of 
Figure [8] (9 w 18°) with the KHs obtained from our F/G- 
dwarf thick-disc subsample, shown in Figures [2] and [5] This 
result clearly indicates that some of the features observed in 
the KH of the Milky Way thick disc, especially the peak at 
E ~ —0.97 x 10 5 km 2 s -2 , could be associated with density 
waves excited by the Sagittarius dwarf galaxy. An important 
implication of the strong azimuthal dependance observed in 
our simulations is that the orbital properties of the Sagit- 
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Figure 7. Top panel: Distribution of particles in the X - Y 
plane, as obtained from the Heavy Sgr simulation after approxi- 
mately 2.1 Gyr of evolution. The different colours (contours) in- 
dicate different numbers of particles. As a consequence of the 
non-axisymmetrical energy kick imparted by the satellite, spiral 
features are induced in the galactic disc. The location of the Solar 
neighbourhood-like volumes analysed in Section l6.2l are indicated 
with solar symbols. The cyan dot indicates the local volume that 
best resembles the observations in the Solar Neighbourhood. The 
white dots labelled as "1" and "2" trace the regions where the 
Heavy Sgr crosses the plane of the disc for the first and second 
time, respectively. Its current projection on the X - Y plane is 
indicated with the symbol "C". Bottom panel: Evolution of the 
galactocentric distance of both Sgr models as a function of time. 

tarius dwarf galaxy could be further constrained by fitting 
the locations and amplitude of the significant peaks in the 
models to those observed in stellar samples. 

In the middle panels of Figure [5] we show the results 
of the previous analysis obtained from the Light Sgr model. 
The Solar neighbourhood-like spheres were placed at the 
same locations. As before, the black solid lines show the 
kernel histograms obtained after computing particle energies 
using the analytic potential, whereas for the red dashed lines 
the A r -body potential was used. As expected from a less- 
massive satellite (and thus a weaker perturber), the peaks 
previously observed in the kernel histograms of the Heavy 
Sgr model are either very mild or not present at all. It is 
interesting to note that, as previously discussed, the Light 
Sgr model results in a much stronger bar than the Heavy 



Sgr model. The more massive satellite suppresses bar for- 
mation as a result of enhanced central disc heating. This 
clearly demonstrates that none of the observed features in 
these simulations could be associated to resonances with 
the galactic bar. Otherwise, we would expect to observe 
higher-amplitude features in the Light Sgr model's kernel 
histograms as a result of the interaction with the stronger 
bar. Thus, models of the Sagittarius dwarf galaxy progenitor 
with a total mass 1O 1O ' 5 M0 could be ruled out. Moreover, 
note that its total mass could be further constrained by com- 
paring the amplitudes of the significant peaks observed in 
the kernel histograms of models with those of the Milky Way 
disc. 

In the bottom panels of Figure [8] we explore Solar 
neighbourhood-like spheres of 1 kpc radius centred at dif- 
ferent galactocentric distances. The spheres are separated 
by 0.5 kpc, and are all located at 6 = 0°. The black solid 
lines show the kernel histogram obtained from the Heavy 
Sgr model. A strong dependence on the amplitude of each 
peak with galactocentric radius can be observed. As we move 
towards the galactic centre we better sample inner den- 
sity waves, represented by peaks at lower energies. Again, 
note the similarities between the KHs obtained in the inner 
spheres (e.g., R — 7.5 kpc) and those obtained from our 
F/G-dwarf sample (see e.g. Fig. [5]). The blue dashed lines in 
the bottom panels of Figure [8] show the KH obtained from 
the Light Sgr model. As expected, mild signatures of density 
waves can be observed only in the outer spheres, where the 
perturbation from the satellite galaxy is stronger. In Figure[9] 
we show the time evolution of the kernel histogram of par- 
ticles located within a given Solar neighborhood-like sphere 
for ~ 2.5 Gyr. To track this volume over time, our sphere 
is rotating with an angular frequency set by the velocity of 
the LSR, obtained after correcting the final mean rotational 
velocity profile for axysimmetric drift. The KH are shown at 
four different times: the initial conditions, the two pericen- 
tre passages, and the present day configuration. From this 
figure it is possible to appreciate the smooth distribution of 
particles in energy space at t = 0. In addition, it is possible 
to observe that strong perturbations in the disc particles en- 
ergy distribution are only excited after the first pericentre 
passage of the Heavy Sgr satellite, as shown by the black 
solid line. Note that the smaller number of peaks in the KH 
of the Heavy Sgr model with respect to the KH shown in 
Figure [T] is due to the shorter period of time this simulation 
has been evolved. The blue dashed lines show KHs obtained 
from the light Sgr model. Notice the approximately smooth 
energy distribution at all times. 



7 SUMMARY AND CONCLUSIONS 

We have analysed the SEGUE F/G-dwarf and the Schuster 
et al. (2006) stellar samples in search for imprints of minor 
merger events in the Milky Way disc. In contrast to a number 
of previous studies, we look for perturbations in the phase- 
space of the host disc rather than merger remnants. We ap- 
ply the same analyses to high-resolution N-body simulations 
of the interaction between the Sagittarius dwarf galaxy and 
the Milky Way. We summarise our results as follows: 

• The alpha-enhanced, likely thick-disc F/G-dwarf sub- 
sample exhibits significant peaks in the stellar energy kernel 
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Figure 9. Kernel histograms of energies obtained from both Sgr simulations at different times. Each panel shows a KH obtained from 
a distribution of particles within a Solar neighbourhood-like sphere located at R = 8 kpc and rotating with an angular frequency set set 
by the velocity of the LSR, after correcting for axysimmetric drift. The KH are shown at four different times: the initial conditions, the 
two pericentre passages and the present day configuration. The black solid lines show KHs obtained Heavy Sgr model. Note the smooth 
distribution of particles in energy space at t = 0. Strong perturbation in the disc particles energy distribution are only excited after the 
first pericentre passage of the satellite. The blue dashed line show the KH obtained from the light Sgr model. Note the approximately 
smooth energy distribution at all times. 



histograms (KH). Conversely, the thin-disc subsample shows 
a much smoother distribution. This is consistent with the 
expectation that merger-induced perturbations in the host 
galaxy disc (e.g., the Milky Way) create structure in the 
stellar phase-space density, which survives until later times 
only in the old, hotter stellar population. 

• We show that, to unambiguously identify merger- 
induced waves (ringing), samples of stars with velocity er- 
rors Var < 11 km s _1 are required. This depends on the 
time since the waves were excited, as well as on the mass of 
the perturber. 

• Complementary evidence of ringing can be obtained 
from the SCH06 catalog. The energy kernel histogram of 
stars here also reveals peaks, most of them remarkably con- 
sistent with those observed in the F/G-dwarf sample. 

• Simulations of the interaction between the Milky Way 
and the Sagittarius dwarf galaxy show that a relatively- 
massive progenitor (i.e., masses of » 10% of the host) can 
induce density waves in the Galactic disc. Furthermore, the 
features seen in the energy kernel histogram of the observa- 
tional data can be qualitatively well reproduced with these 
simulations. This is interesting since, by matching features 
in the kernel histograms of both models and observations, 
one could attempt to further constrain the orbital proper- 
ties of Sagittarius, thanks to the strong azimuthal and ra- 
dial dependance of the distribution of peaks. Note, however, 
that other possible origins for these perturbations need to 
be ruled out before performing this kind of analysis. 

• On the contrary, low-mass models of Sagittarius barely 
excite density waves in the disc. Thus, if we assume that the 
observed peaks in the kernel histogram of the F/G-dwarf 
sample are the response of the Galactic disc to the tidal 
interaction with Sgr, we could rule out such low-mass mod- 
els. Furthermore, by matching the amplitudes of the peaks, 
one could attempt to further constrain the total mass of the 
Sagittarius progenitor. 

Comparing Figures [5] and [8] we note that variations 
in the azimuthal and radial positions of the Sun with re- 
spect to the Sagittarius orbital plane can result in similar 



structure in the stellar energy distribution. This degener- 
acy can be broken by studying in detail the variation of 
features in the stellar phase-space distribution when posi- 
tion in the disc is varied. This in turn requires larger and 
deeper stellar samples, which will become a vailable in ongo- 
ing a nd future su rveys, such as APOGEE (j Maicwsk i et al.l 
l2010i ). HERM ES dWvlie-de Boer fc FreemaiJl201Ch . anduT 
timately Gaia l)Perrvman et al.ll200ll ). 

Although we cannot fully discard the possibility that 
the observed features have some secular origin, the much 
smoother energy distribution of the thin-disc population 
make this scenario unlikely. If these features had originated 
by resonant interaction with, e.g., the Galactic bar, we would 
expect them to be stronger in the thin disc, rather than in 
the thick disc. This result is consistent with the hypothesis 
that features associated with ringing should last longer in 
the thick disc. Thick-disc stars spend relatively little time 
near the Galactic plane, where heating from spiral arms or 
the bar and scattering by giant molecular clouds is most vig- 
orous. In addition, features generated by the bar or spiral 
arms can be observed as lumps in the u-v plane. We have 
explored this plane, and found no significant overdensities 
that could be linked to this kind of secular perturbations. 

In this study we considered simulations with two dif- 
ferent masses for the Sagittarius progenitor, representing a 
lower (Mvit = 10 10 5 M ) and upper (M vir = 10 11 M ) 
limits. We have found that the Light Sgr cannot explain 
the structure seen in the energy distribution of our ob- 
served samples. On the other hand, the Heavy Sgr model 
creates strong disturbances in the Milky Ways disc, result- 
ing in prominent peaks in the energy kernel histogram in 
localised spatial volumes, perhaps even stronger than what 
is seen in the observations (compare Figures [5] and [8} . As 
shown by G12, satellites with masses 10% and 20% of the 
mass of their host generate density waves, which can be fol- 
lowed for no longer than 1 and 2 Gyr after the satellite's full 
disruption, respectively. Close encounters between massive 
substructures and galactic discs should be common occur- 
rences in a ACDM cosmology since z = 1 (Kazantzidis et 
al. 2009). However, the short relaxation time scale involved 
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and the lack of observational evidence of a merger with such 
characteristics in the last 1-2 Gyr suggest that Sgr is the 
most likely perturber. 

It is important to note that all simulations analysed 
in this work considered initially stable, unperturbed stellar 
discs. It remains to be explored how a previously perturbed 
disc, either by a merger event or other dynamical process, 
would react to the tidal interaction with a satellite galaxy 
such as Sgr. Future studies should consider simulations span- 
ning a range of masses for the Sgr progenitor, between the 
two limits we presented here, as well as varying Galactic disc 
models, which can be tuned more finely as our knowledge 
of the dynamics of the Milky Way expands. In addition, it 
remains to be studied whether perturbations from multiple 
satellite galaxies can be successfully disentangled. 
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